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Abstract. We show that the response to an electric field, in models of the Integral 
Quantum Hall effect, is analytic in the field and has isolated essential singularity 
at zero field. We also study the breakdown of Chern numbers associated with the 
response of Floquet states. We argue, and give evidence, that the breakdown of Chern 
numbers in Floquet states is a discontinuous transition at zero field. This follows from 
an observation, of independent interest, that Chern numbers for finite dimensional 
Floquet operators are generically zero. These results rule out the possibility that the 
breakdown of the Hall conductance is a phase transition at finite fields for a large class 
of models. 



1. Introduction 

The principal motivation for the present work is the question: Is the breakdown of the 
Integer Quantum Hall effect a (quantum) phase transition? Since the Hall conductance 
in the adiabatic limit is identified with a Chern number, the question can also be phrased 
as: Is the breakdown of Chern numbers a phase transition? 

Experimentally, || g |, [TJ, p|, g ||, gj| the breakdown of the Hall effect at finite 



driving currents is signaled by the onset of dissipation, and is accompanied by hysteresis 
and complex dynamical behavior. The critical current and voltage depend, in general, 
on the geometry of the system, the temperature and on the magnetic field. It has 
been suggested that a phase diagram for the breakdown resembles the phase diagram 
of superfluidity . 



Is there a theoretical basis for identifying breakdown with a phase transition? 
Naively, one can argue both ways. In a class of models of the Integer Hall effect the 
Hall conductance (at zero temperature and in the limit of linear response) is related to 
a Chern number pl| . Chern numbers, being integers, depend discontinuously, if at all, 
on parameters in the Hamiltonian. So, if the strength of the external electric field was 
just like any other parameters in the Hamiltonian, one would expect the breakdown of 
the Hall conductance to be discontinuous and at non-zero field. This would say that 
the breakdown of the Hall effect is indeed a quantum phase transition at finite fields. 

f E-mail addresses: avron@physics.technion.ac.il and konsz@physics.technion.ac.il 
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On closer inspection one realizes that this line of reasoning can not be quite right. 
For, if it was, then the Hall conductance would remain precisely quantized also for small 
but finite values of the electric field. This would imply that there are no corrections 
to the quantization - not even exponentially small correction. Common wisdom is 
that while there are no power corrections to the integral Hall conductance, there are 
exponentially small corrections JT3J. As we shall explain in detail, and this is going to 
be a key point of our analysis, the strength of the driving electric field (or the driving 
emf) is a special parameter which affects the Hamiltonian and the evolution in a way 
that is structurally different from say, the strength of the magnetic field or the disorder 
potential. 

A simple and common argument, with some experimental support, says that the 
breakdown occurs at finite driving fields so that the critical field, S c , scales like B 3 ^ 2 , 
where B is the strength of the magnetic field. This estimate follows from comparison 
of the energy gap in Landau levels with the voltage drop on a magnetic length. The 
breakdown is then attributed to tunneling between Landau levels. This argument does 
not directly address the question if the breakdown is a phase transition. It also has a 
weakness in that the Landau Hamiltonian with constant electric and magnetic fields is 
explicitly soluble and does not show breakdown. For other theories of the breakdown 
see e.g. [p| and references therein. 

To address the breakdown as a quantum phase transition, a handle on the 
conductance at finite fields is needed. This goes beyond linear response and Kubo's 
formulation. 

In this work we shall concentrate on the breakdown of the Hall conductance, rather 
than the breakdown that occurs in the dissipative conductance in the Hall effect. This 
is done for two reasons. The first is that we are interested in breakdown that occurs in 
Chern numbers, for which the Hall conductance is a basic paradigm. The second is for 
concreteness sake. Some parts of what we say can be transcribed, mutatis mutandis, for 
the dissipative conductance. 

We find no theoretical support to the hypothesis that the breakdown in the Hall 
effect is a phase transition at finite fields. Rather, we find support to the claim that the 
conductance has an essential singularity at zero fields, which can, of course, manifest 
itself in what resembles a phase transition. 

In section ^| we show that the strength of the driving electric field (or emf), S, 
can be related to a time scale r. In section |3| we show that the expectation value of 
a (bounded) observable, in particular, the current density in tight-binding models of 
non-interacting electrons, is an analytic function of S with an essential singularity at 
£ = when St is kept fixed. In section ^ we consider the analytic properties of a 
natural notion of transport for a class of model Hamiltonians which, in the limit of 
linear response, reduces to the usual notion of conductance that coincides with a Chern 
number. We show that this observable is analytic in S with an essential singularity 
at S = 0. It is interesting that the absence of phase transition at finite fields can be 
shown for the same class of models where one can prove quantization. In section [5| we 
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study the Harper model for which we present numerical results. In section [7| we study 
Chern numbers associated with Floquet states, their properties and their interpretation 
as Hall conductance of the Floquet states. In section |8| we study the breakdown of Chern 
numbers associated to Floquet states. We show that the breakdown is discontinuous and 
argue that it occurs at zero fields, £ = 0. In section |9] we describe numerical evidence 
that supports the claim that non-zero Chern numbers for Floquet states are unstable 
against perturbations in the Hamiltonian. 



2. Driving Fields Interpreted as a Time Scale 

Because gauge invariance allows one to impose one condition on the scalar and vector 
potential, it is always possible to choose a gauge where the external electric and magnetic 
fields are described only by the vector potential A. In order to produce electric field 
this potential has to be time dependent. Assume that this dependence is characterized 
by some time scale r as A(x,t/r). Suppose now that we scale time so that t = ts. 
Maxwell equation gives the external electric field, (in scaled time), as 

E(x,s) = -^d s A(x,s). (1) 

With A fixed, weak electric fields correspond to large r while strong electric fields 
correspond to small r. In systems that are otherwise time independent, r interpolates 
between weak and strong fields. A similar argument can be made about the emf, which 
is a line integral of the electric field. 

The identification of the time scale r with the strength of the external driving, be 
it the electric field or an emf, is not a new idea, of course. It lies at the heart of the 
identification of the adiabatic limit with linear response. What is perhaps new here 
is that we want to use this correspondence for any driving, and in particular identify 
strong driving fields and large emfs with short time scales. Equation ([!]) suggest that 
we write 

where £ is a measure of the strength of the driving field. We have put fundamental 
constants into this relation so that (in cgs units) £ has the dimensions of an emf (or 
voltage). Because lattice models also have a natural length scale - the lattice spacing, 
one can choose constants so that £ has dimensions of an electric field. The identification 
of the strength of the field with an inverse time scale is central to our considerations 
and turns out to have consequences for transport. The first and easy consequence is 
that the question of phase transition in, say, the Hall conductance as function of the 
driving electric field, can be phrased as a question about the analytic properties of the 
Hall conductance as a function of r. The second consequence is discussed in the next 
section. 
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3. Analyticity of Observables in r 

The identification of the time scale with the strength of the driving makes the driving 
field, be it an electric filed or an emf, a special parameter in the Hamiltonian. This can 
be seen from the form of the Schrodinger equation for the time evolution operator in 
scaled time: 

iU T (s) = TH(A(s))U T (s), U T (0) = 1. (3) 

The r dependence of the Hamiltonian is linear, and a-fortiori, analytic in r, irrespective 
of how H depends on A. This has the consequences: 

Proposition 1 Suppose that H(A(s)) is bounded and self- adjoint, thenU T (s) andUl(s) 
are both entire functions ofr. 

Proof. The first part follows from a standard argument about the absolute convergence 
of the Dyson series for U T . The second assertion is a consequence of the fact that, for 
real r, UJ. satisfies 

iUl(s) = -rUl(s)H(A(s)), flj(0) = 1 . (4) 

The Dyson series implies that one can extend UJ.(s) to an entire function of r. ■ 
It follows that: 

Corollary 2 Suppose that I is a (fixed) bounded operator, and ip a (r independent) 
initial state then: 

I(T,s) = (il>\l4(s)IU T (s)\l>) (5) 
is an entire function of r. 
This leads to the main result of this section: 

Theorem 3 Let H(A) be bounded self-adjoint operator. Let I be a bounded observable, 
then, its expectation value (I)(£,s) is an analytic function of £ with isolated essential 
singularity at £ = and with Laurent expansion 

W(*>*) = £^> ( 6 ) 

n=0 

with infinitely many of the a n ^ (except if (I) is a constant). 

Proof. From proposition III, (I)(£,s) is an analytic function whose Taylor expansion 
about r = is absolutely convergent with radius of convergence that is infinitely large. 
The expansion can not have positive powers of £ , for if it did, the response at r = 
would not be analytic. The singularity at the £ = must be essential for if was a pole 
then the response would diverge as r — ► oo along any direction. But, the response is 
bounded on the real r axis by the self-adjointness of the Hamiltonian. It follows that 
the response can not have a pole of finite order at £ = 0. ■ 
Remarks: 
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(i) It is important for the conclusion of the theorem to hold that the scaled time s is 
kept finite. If one lets s —>■ oo then, analyticity in r may be lost. 

(ii) We have restricted ourselves to bounded Hamiltonians and bounded observables. 
This is because phase transitions are normally a long wavelength, low energy 
phenomenon. Some parts of the theorem can also be extended to unbounded 
Schrodinger type Hamiltonians provided some care is taken about questions of 
domains of operators. 

(iii) The smooth dependence of the evolution in r is not really special to quantum 
mechanics. It holds also in classical mechanics under slightly stronger conditions, 
namely that d p H and d x H are both bounded functions (provided the initial state 
is analytic in r). 

(iv) A classical model that shows breakdown of analyticity in a constant electric field is 
the washboard potential with initial state at rest at a local minimum. When £ is 
sufficiently small, the velocity stays zero for all times. When £ passes a threshold 
the particle accelerates indefinitely. This is not a counterexample to the analyticity 
of observabales because the initial state in the washboard potential is £ dependent 
in a non- analytic way. 

(v) Two prototype functions for (/) that satisfy the conclusion of the theorem are 
exp(— c/£ 2 ) and sin(c/£ 2 ). 

(vi) The theorem has an interesting implication to the question of power law corrections 
to linear response. It follows from the theorem that if (I) has an asymptotic 
expansion in powers of £ G IR+ at £ = then this expansion must vanish identically. 
One may wonder if the absence of power law corrections to linear response is a 
valid conclusion of equation (||). It is not, as one can see from the second prototype 
example in (v), where an asymptotic expansion in powers of £ G M + does not exist. 

(vii) Absence of power corrections to linear response in quantized Hall conductance has 
been proven rigorously for the quantum Hall effect by Klein and Seiler |13 |. Their 
proof uses an adiabatic theorem to all orders, to show that an asymptotic expansion 
exists and a clever trick that shows that the coefficients must vanish. The theorem 
above can be used to replace their clever trick. 

In tight-binding models the current operator, for finitely many interacting electrons 
or the current density operator for infinitely many non-interacting electrons, is a 
bounded. It follows that the currents in such models have an essential singularity at 
£ = 0, but are analytic at non-zero £ provided the time t£ is kept fixed. 



4. Breakdown of Chern Numbers 



In this section we consider the breakdown of an observable associated with quantized 
charge transport. Quantized charge transport occurs for a class of Hamiltonians in the 
adiabatic limit. This class of Hamiltonians includes models of the quantum Hall effect, 
and in particular includes the Harper model. It also includes certain models of the Hall 
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effect with electron-electron interactions. The observable that we consider reduces to 
the Hall conductance in the limit of linear response, and coincides with a Chern number. 
In this section we discuss the analytic properties of this observable with r. 

The model Hamiltonians for which quantization occurs have the following struc- 



ture |2(], H((f), k) is a self-adjoint Hamiltonian that depend periodically on two real, 
dimensionless parameters, and k, with period 2tt. if(0, k) may be associated with 
a finite multiparticle system, where and k are external parameters, for example, two 
Aharonov-Bohm fluxes. Alternatively, H((j), k) may be a Bloch type Hamiltonian in two 
dimensions describing infinitely many non-interacting electrons, where and k are two 
Bloch momenta. In either case, we shall require that for fixed and k the Hamiltonian 
H((p, k) has discrete spectrum with no eigenvalue crossing. For the sake of simplicity we 
assume that the and k dependence is smooth and that H (0, k) is a bounded operator 
such as a tight binding model and its multiparticle generalizations. 

The time dependence comes from a time dependence of on a time scale r[|. We 
suppose that 0(s) is a smooth, monotonically non- decreasing function of s with <ft(s) = 
in the past, s < 0, and 0(s) = 2ir, in the future, s > 2tt. H is therefore time dependent 
only on a finite interval of (scaled) time [0, 2n]. 

Since d^H is the current operator in these models the total charge transported by 
the action of is: 

Q(r; V) = ^ jT dk As (u T (s, k)ij(k) dH{ ^ ] ' k) U T (s, k)m) , (7) 

with initial condition ip that is an eigenstate of H(<p, k) for s = 0. 

There are several special things that happen in the adiabatic limit. First, Q 
coincides with the Hall conductance defined via Kubo's formula |2T| . Second, it 
is independent of the functional form of (provided satisfies the limiting conditions). 
Third, since in the adiabatic limit there is no current once the driving stops, i.e. when 
= 0, Q can also be written as 

provided s > 2n. Q is always a measure of the charge transport, but its identification 
with a conductance is valid, in general, only in the adiabatic limit|[ 

Applying the results of the previous section we see that the breakdown of the Chern 
number Q is smooth, and has no phase transition at finite fields. Due to the prefactor 
r in equation @, ao = in equation (|6|). 

It is interesting that for the class of models where one can prove that Q is quantized 
in the adiabatic limit, one can also show that it is an analytic function of the field away 

X Here 4> plays the role of A of the previous section. 

§ Equation (|^) implies that Q vanishes in the limit r — > 0, i.e. in the limit of large external fields, 
contrary to common experience. One reason for this is that tight binding models are unreasonable 
when the external fields are large on atomic scale. We shall take the point of view that the breakdown 
is a low field phenomenon that is divorced from the asymptotic behavior at very large fields. 



Q( S ;r,^) = ^^ dA; J ds' (u T if> 
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from £ = 0. (For model Hamiltonians with infinitely many interacting electrons there 
is, at present, no proof of quantization either.) 

There are now three possibilities. The first, and perhaps simplest, is that the 
breakdown of the Hall effect is a consequence of an essential singularity at zero field. 
The second is that the breakdown of the Hall effect is associated with the limit s —>■ oo. 
And the third is that the breakdown is a property of infinitely many interacting electrons. 
We shall examine the second possibility in section [7|. 



5. Example: Hall Conductance in the Harper Model 

The Harper model is the simplest, non-trivial, model where one can study the breakdown 
of the Hall effect in detail, at least numerically. The model is associated with a square 
lattice, Z 2 ; an external homogeneous magnetic field B, and homogeneous electric field £ 
pointing in the x direction. We choose a gauge so that the electric field is described by 
a time dependent vector potential. After separation of variable, the model is described 
by a Hamiltonian on Z parameterized by one Bloch momentum, k. The Hamiltonian 
action on the vector \1/ G ^ 2 (Z), while the electric field is acting, is: 

e i£t y x+1 + e~ i£t ty x „ l + 2cos{Bx + k)^ x , fce[-7r,7r], x G Z . (9) 

For rational magnetic field B = 2-npjq with p,g 6 Z the Hamiltonian is periodic in 
x, with period q. One then classifies the solutions by a second Bloch momentum, 
£ G [— 7r, 7r] so that ^/ x+q = exp(— i£)^f x . Fixing periodic boundary conditions is achieved 
by the unitary transformation ty x — > e^ x ^/ x . So, finally, the requisite form of the Harper 
Hamiltonian we shall study is: 

(H(<j>, fc)*) = e^+x + e-***.! + 2 cos(2tt^ + k)* x , 



m x+q = V x , <j>{t) 



q 

if t < 0; 



q 



£t + -, if0<t<#; (10) 
q c 

£ 

2tt H — , otherwise. 

q 

This corresponds to a qxq hermitian matrix, periodic in and k. The driving electric 
field, £, is related to the adiabaticity parameter r, by (§) (in units where e = h — 1). 
In this example 0(s) is continuous and piecewise linear. As a consequence the electric 
field is discontinuous in time. It is easy to modify the model so that the electric field is 
switched on and off continuously. We have examined also such models and the behavior 
is qualitatively similar to models with discontinuous switching. 

For p = 1, q = 3, the Harper Hamiltonian is a 3x3 matrix. Its Chern numbers are 
{3, 3, —6}. Figure [I] shows the charge transport, Q(2tt; t, ip), as function of the applied 
field. The graph has a rich and complex structure, but no sharp breaking. Substantial 
deviation from integral quantization occur near r = 5. 
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6. Long-time limit and Floquet states 

In this section we consider the long-time limit of observables, in time-periodic (finite 
dimensional) Hamiltonians. An example is the current density operator in tight- 
binding models, such as the Harper model, driven by a time-independent electric field. 
One advantage of using a time dependent representation over the time independent 
representation is that one has to deal with finite dimensional matrices. In the time 
independent representation the matrices are infinite dimensional. 

In the theory of time-periodic Hamiltonians Floquet states |7], 16, [| play a role 



analogous to that of eigenstates for time independent operators. We shall see that 
the long-time (Abelian) limit of observables in time periodic systems is related to the 
expectation value of the observable in Floquet states. 

Consider a time-periodic, self-adjoint, finite dimensional matrix Hamiltonian H(s + 
2tt) = H(s). The Harper model of section |5] is an example except that now the electric 
field is time independent for all (positive) times s > 0. 

Let F denote the unitary evolution over one cycle 

F = U(2ir). (11) 
The time evolution for n G Z periods is clearly 

F n = U(2mx) . (12) 
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We assume that F is a finite dimensional, non-degenerate, matrix. Its spectral 
representation is then 

F = Y. eiEn \€)(€\- (13) 

rt 

E n are the quasienergies, and \if)^) the eigenvectors of F. 

Consider the observable associated to the bounded operator / (e.g the current 
density operator in tight binding models). Let if) be the initial state of the system at 
s = 0, then the Abelian, long-time, limit of the expectation value of J, when if) evolves 
according to the Schrodinger evolution, is: 

M 1 M 

lim — V(I(2tt7')) = lim — (F j ip 1 1 1 F j ifj) 

J=l J=l 

= £ « I 7 I V? ) <W> M"> ( j™, ^ £ ^^"^ ) 

l,n \ j=l J 

= £KM)rWl J M p > • (14) 

The long time behavior is a weighted sum of the expectation value of the operator / in 
Floquet states, i.e. (ipf \ I \ ipf). 

In particular, if one now defines the Hall conductance as the ratio of the Hall current 
density to the electric field in the long time limit, then the breakdown of the Hall effect 
is related to analytic properties of the expectation value of the current density operator 
in Floquet states. This is the subject of the following sections. 



7. Transport in Floquet States 



In this section we recall, and extend, results of Ferrari [ID that relate the Hall 
conductance of Floquet states to their Chern number and the winding numbers of their 
quasienergies. 

Consider a time-periodic, self-adjoint, finite dimensional matrix Hamiltonian H(s + 
27r, k) = H(s, k), which depends analytically on s and k. The Harper model of section 
[| is an example except that now the electric field is time independent for all (positive) 
times, s is, as before, the scaled time. 

We define the Floquet operator as the unitary evolution over one cycle 

F T {s,k) = U T {s + 2ti, k)Ul{s,k) . (15) 

Floquet theorem can be expressed as 

F T (s + 27r,k) = F T (s,k) . (16) 

We assume that F T is a finite dimensional matrix. It therefore has discrete spectrum 
and its spectral representation is 

F T {s, fc) = E e iE " {s ' k;T) V n (s, k- r) . (17) 
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E n are the quasienergies, and V are the eigenprojections. We shall denote by \ip^) a 
unit eigenvector of F T . 

Lemma 4 For a Floquet operator generated by a bounded and analytic H(s, k), 

(1) The quasienergies are independent of s, E n (s,k;r) = E n {k\r). 

(2) The eigenf unctions (eigenprojections) obey the Schrddinger (Heisenberg) equations 

id a \ifi(s,k)) =TH( 8 ,kM(s,k)) , (18) 
id s V(s, k; r) = t[H(s, k),V(s, k; r)} . (19) 

(3) F T (s, k) is an entire function of r and analytic function of s and k. 

Proof. The Floquet operator satisfies 

F T (s, k) = U T (s, k)F T (0, k)Ul(s, k) , (20) 

so its eigenvalues are independent of s. The eigenfunctions and eigenprojections 
transform in the same way 

\tf(s,k)) = U T (s,kM(0,k)), (21) 
P(s, k; t) = U T (s, k)V(0, k; r)Ul(s, k) , (22) 

from which the second assertion follows. The third item follows from proposition |I|. ■ 
F depends on s, r and k. Because of equations ( p0|) the s variable is uninteresting. It 
is a basic fact of perturbation theory ||rTJ that if a (normal) operator depends analytically 
on two variables (or more), one can choose, in general, at most, one variable so that 
the eigenvectors and eigenvalues are real analytic. We choose real analyticity in the k 
variable. The price one has to pay for this is two fold: First, the r dependence of the 
projections will, in general, be discontinuous at crossings^] and the second is that the 2ti 
periodicity in k may be lost. (One can not assume that V is both real analytic and 2ir 
periodic in k.) But, it is 2nN periodic for some finite N. If there are no crossing then 
N — 1. 



A result of Ferrari |T0[ identifies the Chern numbers of Floquet eigenstates with the 
charge they transport and relates the Chern number of Floquet states to the winding 
number of the quasienergies. This identification will play a key role in the next section 
where we discuss breakdown. 

Proposition 5 Let F T (0,k) be a finite dimensional Floquet operator, which is unitary 
for real values of r and k, and analytic in both. Fix r and label the spectrum so that 
V T (0, k) is real analytic in k with period 2ttN . Then 

(1) The quasienergies E(k;r) and the associated eigenfunctions \ip^(s,k)) are simul- 
taneously real analytic in k and s. 

|| This can be seen already for the 2x2 matrix function ra x + ka z where a are the Pauli matrices. 
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(2) The conductance associated to non-interacting electrons that fill a k-band of Floquet 
states is a Chern number multiple ofl/N: 



2ttN 



Im 



2- 



ds 



2ttN 



dk{ 



ds dk 



(23) 



(24) 



(3) Q(T;ip T ) is also the winding number of the corresponding eigenvalue, 

(4) Q{r,^) = for all£ > 2 max Site[0> ] \\H(s,k)\\. 

Proof. The first assertion follows from the real analyticity of the projection, (|2T| ) 
and standard facts from perturbation theory [TL\. The second assertion follows from: 



T 



■ij'r 



dH 
~dk 



d 

T—(i; T \H\i} T ) + Im( 



dk ds 



(25) 



The observable on the left hand side of this identity is the ratio of the current to 
the driving electric field (in each k channel). This means Q is a conductance. The 
identity fl25|) is a consequence of integration by parts, the equation of motion, and the 
periodicity in k. The quantization of Q follows by integrating over the period in k and 
using standard facts about Chern numbers. Using the identity 



f! 



dF T 



T dk 



(0,k) 



-IT 



2tt 



U}(a,k) 



(dH 



\ dk 



U T (s, k)ds 



(26) 



and equation (Jl 



V? (0, k) | (F%F T ^ (0, k) | V£(0, k)) = id k E(k; r) (27) 

gives assertion 3. Item 4 follows from the observation that the winding vanishes if 
max \E(k; r)| < tt. The maximum can achieved when all the terms in the Dyson 
expansion of the evolution operator have the same eigenstate for their maximal 
eigenvalue and all the eigenvalues are summed up with the right signs. It follows that 

\E\ < 27rrmax||#|| , (28) 

which implies item 4. ■ 
Finally, we note that for the adiabatic limit, £ = 0, eigenstates of the Hamiltonian 
are also Floquet eigenstates. This shows that for £ = the Chern number for Floquet 
states is the Hall conductance of linear response when H(s, k) is the Hamiltonian for 
the Hall effect. 



8. Breakdown in Floquet States 

In the previous section it was shown that the Hall conductance in Floquet states, for any 
field strength £, is related to the Chern number of the state. Since the Chern numbers 
are integers, this implies that the breakdown in Floquet states is always discontinuous: 
a quantum phase transition. For £ = the conductance for Floquet states coincide 
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with the Hall conductance of linear response but when 6 is large enough Q vanish. This 
forces a discontinuous breakdown. 

At first, it appears that one can argue that a breakdown must occur for some finite 
value of 8 because Chern numbers do not change under small deformations of the bundle 
of eigenstates, and so if the Chern number is non-zero at zero field, should it not be 
non-zero also for small fields? 

To understand how the breakdown occurs we first make the observation that non- 
zero Chern numbers for Floquet operators always come with level crossings: 

Theorem 6 Let F T (0,k) be a finite dimensional Floquet operator, which is unitary for 
r, k G R ; and analytic in both. Label the spectrum of F T so that V T (0, k) is real analytic 
in k with period 2tcN. If any Chern number is non-zero then F T has eigenvalue crossing 
for some value of k. 

Proof. By general principles, the sum of all the Chern numbers for any finite 
dimensional matrix, such as F T , must vanish. Hence, if there is a positive Chern number 
for one of the states of F T , there must also be a negative Chern number for some other 
state. The quasienergy with a positive winding must then cross the one with negative 
winding. ■ 
This leads to a puzzle that has Chern on one side, and Wigner von-Neumann 
(WVN) on the other: Wigner and von-Neumann [|I7|, [TJj] say that eigenvalue crossing 



tend to be unstable, unless the dimension of parameter space is three or more^J. The 
relevant parameter space for the Floquet operators is the r-k space which has dimension 
two. This implies that crossings are unstable. On the other hand Chern numbers give a 
topological characterization of the bundle of eigenstates and are stable under continuous 
deformations of the bundle. Who wins? 

The winners, we claim, are Wigner and von-Neumann. It is, of course, true that 
continuous deformation of the bundle keep the Chern number fixed. But, there is no 
reason why a continuous deformations of the Hamiltonian or a variation in r, should 
lead to continuous deformations of the bundle of eigenstates. By perturbation theory 
a small deformation in the Hamiltonian can result in a discontinuous change of the 
bundle at points of crossing. Since non-zero Chern number for Floquet states always 
comes with eigenvalue crossing, the bundles always lie at the boundary of the region of 
stability. As a consequence, generically, at least, Chern numbers for Floquet operators^ 
should be zero and the breakdown should therefore occur at £ = 0+. 

Experimentally, the breakdown of the Hall effect is continuous and consequently, 
the breakdown of Floquet states does not appear to be an appropriate theory for the 
breakdown of the Hall effect. This suggests that the breakdown is not a large time 
phenomenon. Nevertheless, the breakdown in Floquet states is an interesting chapter 
in the theory of breakdown of Chern numbers. 

In the complex case, which is the case relevant here. 
+ At least, those of the class we study here, that depend on a single variable k 
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Figure 2. A map from the 2-torus (ip, 9) to the sphere (the unit vector pointing from 
9 to ip). The linking number is the degree of the map. 



9. Unstable Chern numbers 

We have seen in theorem |6| that non-zero Chern numbers for Floquet operators always 
come with eigenvalue crossings. The question whether Floquet operators do or do not 
have stable non-zero Chern numbers, is equivalent to the question whether Floquet 
operators do or do not have stable eigenvalue winding and crossing. 

9.1. Fragile Winding 

A simple example, that we owe to B. Simon, illustrates how topological objects that are 
normally associated with stability become fragile in the context of eigenvalue problems. 
Consider the unitary 



F has unit determinant so its two eigenvalues are determined by a single angle 0. Clearly 
cos0 = cos(l/r) cosfc. For 1/r = im the two eigenvalues have winding numbers ±1 as 
k goes through a period. But, for most value of r the two eigenvalues repel at and 
7r, and have zero winding. Winding numbers that arise from an eigenvalue problem are 
not stable. 

9.2. The Paradigm 

The basic paradigm of line bundles associated with Chern numbers ±1 comes from 
Berry example of spin 1/2 in a magnetic field. Berry's spin 1/2 model describes a 2x2 
Hamiltonian that is parameterized by the two sphere. In the study of Floquet operators 
and adiabatic transport, one is interested in families parameterized by the 2-torus. The 
question is how to pick a family of 2x2 Hamiltonians on a 2-torus with Chern numbers 
±1. A nice geometric way (figure |2|) to think about the family is to consider the pull- 
back of Berry's spin 1/2 on the 2-sphere, by the Gauss map from the 2-torus to the 
2-sphere with degree oneQ. 

* We thank I. Klich for pointing this out. 




(29) 
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A useful fact about the quasienergies of 2 x 2 Floquet operators with determinant 
1 is that Ei = —E 2 . Hence, gap closure occur only at E — and E = ±7r. A stable 
non-zero Chern number then requires a stable crossing at both at E = and E = ±7r. 
Since this feature simplifies the numerical analysis considerably the 2x2 Hamiltonian 
on the torus should also be traceless. 

An example of a periodic, traceless 2x2 matrix Hamiltonian with Chern numbers 
±1 is 

*(..*)= ; 1 ( MS) " . ,A (30) 

/ . / s 2 , i a \ a - sin s / 
y sm(sy + \a\ \ v ' J 

with 

a = cos(fc + j) + i cos(fc - -) + (1 + i) (cos(s) + 1) . (31) 

The Floquet operator associated to this Hamiltonian has to be calculated by 
numerical integration of the Schrodinger equation (Q). This Floquet is our basic 
paradigm. Now, if non-zero Chern numbers for Floquet operators are generically 
unstable, this will manifest itself in the opening of gaps at 7r and 0. The gaps are 
shown in figure figure ^ and are consistent with the claim that Chern numbers of the 
Floquet operator are generically zero. 




Figure 3. Gaps width of the Floquet operators from Hi. The solid line is for the gap 
near E = tt, the dashed line is for the gap near E = 
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9.3. Another Toy Models 

We have also carried out extensive numerical studies of the following toy model of a 2x2 
traceless periodic matrix family: 



This Hamiltonian is inspired to some extent by the Harper model on a triangular lattice 
HJ. It has the following features: 

(i) For all s,k,e ^ ±1 the spectrum is constant Spec(H 2 (s, k,e)) = {—1, 1}. 

(ii) The Chern numbers associated with the eigenvectors of H 2 {s,k,e) are ±2 for 



We have added an extra parameter, e so that the associated Floquet operator, F T (0, k; e), 
depends on three variables, {r, k, e}. The reason for including e shall become clear below. 

What should one expect for the Floquet operator based on WVN genericity 
argument? Since (the unitary) F T (0, k;e) depends on three real variables, one expects 
isolated points in {r, k, e} space where there is crossing. In other words, for some special 
values of e one expects to find crossing, while for most values of e, there should be no 
crossing, and all Chern numbers should vanish. This turns out to be the case. 

Here too the Floquet operator has to be calculated by numerical integration of the 
Schrodinger equation (Q). 




(32) 



a = cos(fc) , b = cos(s) 



c = cos(s + k) . 



1 < s < 1 and zero otherwise|} 








2 



4 



6 



8 



10 



Figure 4. Plot of the minimal quasienergy gap vs. r for e = 0. Notice the logarithmic 
y-scale. The calculation accuracy was of order 10~ 8 . The circles represents the 
calculation points which are spaced non-uniformly. 



(t This can be seen by examining sections of the bundle. 
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For e = the crossing at E = is stable for all r. The gap at E — ±7r is plotted 
in figure |j. As can be seen, for most values of r the gap is wider then the numerical 
tolerance. This gap is usually very small. It is easy to mistake these small gaps for 
crossing. 

Near integers values of r the gap width decreases rapidly. At these points the gaps 
seem to close, up to our numerical accuracy, and the winding numbers change from zero 
to ±2. 

These results are in conflict with both the putative stability of Chern numbers, 
and the WVN instability of crossing. It is in conflict with Chern, as non-zero Chern 
numbers appear to occur at isolated points on the r axis rather than intervals. It is also 
in conflict with WVN, for, by the no-crossing rule one expects no crossing in the r — k 
plane at all. The model would be consistent with WVN if the point e = turns out to 
be a special point. Although we do not know if and why e = is special, one expects 
WVN to fail for some e, and a test would be to see if the crossing disappear when we 
wiggle away from e = 0. 

For non-zero values of e both gaps closure at E = 0, n indeed disappear, as is seen 
in figure |5], and Wigner and von-Neumann are vindicate. 




sec 



For more numerical results and description for some of the numerical procedures 
13. 



In conclusion, the numerical results for both toy models are in agreement with 
WVN and with our argument that Chern numbers for Floquet operators are generically 
zero. 
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9.4- The Harper Model 

We have also studied numerically the question of stability of Chern numbers for the 
Floquet operator that is associated with the Harper model. 

Telling a true crossing from an avoided crossing in a numerical study is often a 
challenge. Often one finds complicated pictures where it is difficult to tell with certainty 
which is really the case since even real gaps tend to be very small. This is illustrated 
by numerical results for the Harper model with p = 1 and q = 3,4, see figures § and fj]. 

(b) x=2.5 

-"^ ' — 1 [T/ T! 

A 

■ 

^ 2 , , C ,_ 

2 k 4 6 




Figure 6. The quasienergies E(k) of the Floquet operator associated with the 3x3 
Harper Hamiltonian for r — 1.5, 2.5. 



LU 



3 
2 
1 

-1 
-2 
-3 







(a) r=1.5 



2 k 4 6 



(b) t=5.0 



A V / Y* t / f V / " 



2 k 4 



Figure 7. The quasienergies E(k) of the Floquet operator associated with the 4x4 
Harper Hamiltonian for r = 1.5, 5.0. 



Our numerical results are consistent with the assertion that the Harper model has 
some stable crossing. It is an open problem to prove or disprove this. 
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10. Failure of Berry's method to detect Floquet crossings 

One reliable method to identify crossing of self- adjoint matrices, is to compute a Chern 
number (in the Hermitian case) or the Longuet-Higgins phase (in the symmetric case) 
|2]. |3| . Unfortunately, this method does not work for the kind of Floquet operators that 
we consider, as we now proceed to explain. 

Because of the s independence of the spectrum an isolated crossing would look like 
the line of constant k and r as in figure || The associated Chern number involves the 
surface of integration, a 2-torus of (s,6), as shown in figure || 



▲ 






\ Line of 
degeneracy 






e,,i 


k 








— ► 








Integration 



Figure 8. The 2-torus which is used to enclose the crossing line. Notice that s is a 
cyclic parameter, so the cylinder closes to a torus. 



Using similar arguments as in equations (^3|) and (|24|) we find 

4lP) .4f W)fc -*hM, ,33) 

Since the surface lies away from any crossing the winding number of all eigenvalues 
must be the same. Since the Chern numbers sum to zero, the winding must be zero too. 
Chern number are therefore not useful to identify the crossing of Floquet operators. 



11. Conclusion 



Transport coefficients in linear response theory, such as the Hall conductance, can, in 
general, undergo a phase transition as functions of parameters in the Hamiltonian. This 
is because, in the linear regime, the conductance can be expressed in terms of the spectral 
data of the Hamiltonian. Spectral data may, and in general will, lose smoothness near 
eigenvalue crossing and allow for the loss of analyticity. 

The breakdown in Chern numbers associated with the Hall conductance in finite 
fields seen in figure [l] is not associated with loss of analyticity at finite fields but with 
an essential singularity at zero fields. 
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One way to have Chern numbers undergo a phase transition with the applied field 
is to define them by a spectral problem. This is the case for Chern numbers of Floquet 
states. An external time independent electric field, in an appropriate gauge, leads 
to a tight-binding Hamiltonians that is periodic in time and admits Floquet analysis. 
The Hall conductance in Floquet states reduces to a spectral problem for the Floquet 
operator. This holds for any value of the driving field, and not just in the linear response 
regime. Because of this, changing the field can lead to a non-smooth behavior of the 
conductance when eigenvalues of the Floquet operator cross. We argue that the Chern 
numbers of Floquet operators are generically zero, because non-zero Chern numbers 
come with eigenvalue crossing. Non-zero Chern numbers are unstable to small variations 
in the Hamiltonian and the breakdown occurs at zero fields. 

In conclusion we can say that we find no theoretical reason for the breakdown of 
the Hall effect to be a quantum phase transition at nonzero fields, although we do not 
rule this out for more complicated models such as models of infinitely many interacting 
particles. In addition we observe that under the condition where one can show that the 
charge transport is quantized, one can also show that it is an analytic function of the 
field with an essential singularity at zero field. 
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